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A survey is given on the present knowledge of the polarized parton distribution func- 
tions. We give an outlook on further developments desired both on the theoretical as 
well on the experimental side to complete the understanding of the spin-structure of 
nucleons in the future. 

1 Introduction 

Deeply inelastic scattering provides a clean way to extract the parton densities of nucleons. 
After the initial observation that the nucleon spin is not formed by the quarks dominantly [1], 
detailed measurements of the polarized structure functions followed during the last 20 years. 
The central question concerns now the parton distribution functions and their scale evolution 
rather than just their first moment. Since the nucleon spin receives also contributions from 
the angular momentum of the quarks and gluons, these degrees of freedom have also to be 
studied. This requires the analysis of non-forward scattering cross sections. In inclusive 
deep-inelastic scattering the sensitivity to resolve the different sea-quark contributions is 
rather limited. One way to extract this information consists in measuring semi-inclusive 
processes [2]. A central question concerns the polarized gluon distribution, which can be 
determined from the scaling violations of <?i(x, Q 2 ), deep-inelastic heavy flavor production, 
and hard processes measured at hadron colliders. The inclusive and semi-inclusive hard 
processes in polarized scattering provide an important laboratory to test QCD. More than 
exploring the level of twist-2, which has been investigated in some detail already, one may 
probe the twist-3 contributions in various transverse spin processes. In the following we 
survey the theoretical status of inclusive polarized deeply inelastic scattering and the status 
of polarized parton densities, cf. also [3]. We close with an outlook on investigations required 
in the future. 

2 Theoretical Aspects 



In the deep-inelastic domain the polarized nucleon structure functions receive contributions 
of leading and higher twist, depending on the region in Q 2 and W 2 probed. The leading 
twist contributions are those of twist r = 2 for gi(x, Q 2 ) □ and r = 2, 3 for g2(x, Q 2 )- At 
the level of twist-2 one may extract the polarized parton densities from the data on the 
structure function gi(x, Q 2 ) performing an analysis in the framework of perturbative QCD. 
During the past decades higher orders have been approached steadily. The running of a s (/x 2 ) 
is known to 0(af ) [5], both the polarized anomalous dimensions [6] and massless Wilson 
coefficients [8] were calculated to 0(a 2 ) and the first non-singlet moment, the polarized 
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a Note that gi(x,Q 2 ) contains twist r = 3 contributions due to target mass corrections, [4]. 
b Due to the Ward identity = AP^| this splitting function is also known to O(ctf), [7]. 
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Bjorken sum-rule, to O{o 3 s ) [9]. The heavy flavor Wilson coefficients, in the whole kinematic 
region, are only known to 0(a s ) [10]. For Q 2 ^> m 2 , i.e. in the region Q 2 /m 2 <; 10, the 
Wilson coefficients were calculated in 0(a 2 ), [11]. An interesting property is exhibited by 
the gluonic heavy flavor Wilson coefficient, the first moment of which vanishes in leading 
and next-to-leading order [10,11]. Given a positive polarized gluon density, this implies a 
negative correction to gi(x, Q ) in the region x i5 10 -2 and a positive contribution above. 
Conversely, the Wandzura-Wilczek relation implies a positive correction for the twist-2 
part of g2(x, Q 2 ) for x 2 ■ 10 -2 , but a negative correction for larger cc-values, cf. [12]. The 
anomalous dimensions for the evolution of the transversity distribution are known to 0(a 2 ) 
in general [13] and for a series of moments to O(a^) [14]. At present only next-to-leading 
order QCD analyzes can be performed to extract the polarized parton distributions. 

At the level of the twist 3 contributions to the polarized structure functions several sum 
rules and integral relations were derived, cf. [4,15]. The twist-2 contributions to the structure 
function g2(x, Q 2 ) is given by the Wandzura-Wilczek relation [16] for the quarkonic, gluonic, 
heavy flavor contributions, target mass corrections, and even diffractive scattering, cf. [4,12, 
17]. The twist-3 contributions to the structure function g 2 (x,Q 2 ) were calculated to one- 
loop order. The 0(a s ) non-singlet and singlet anomalous dimension matrices, respectively 
their corresponding expressions in momentum fraction space, were derived in Refs. [18] using 
different techniques. The 0(a s ) Wilson coefficients were calculated in [19]. Although the 
precision on g2{x,Q 2 ) improved during the last years [20] and some difference between the 
data and the Wandzura-Wilczek approximation is seen, still more precise data are required 
to test the QCD-predictions. First non-singlet moments of the twist-3 operator expectation 
values were determined in lattice simulations [21,22]. 

Also in case of the polarized structure functions small- a; resummations are discussed, 
which can be described on the basis of infrared evolution equations [23] and emerge both 
for the non-singlet [24] and singlet structure functions [25,26]. These resummations apply 
to the leading pole 0((a s /iV 3 ) fc )-terms only, with N the Mellin- variable. Performing the 
rcsummation one obtains a branch-cut instead of the perturbative pole-terms, which yields 
a milder singularity. The rcsummation is consistently accounted for by the Callan-Symanzik 
equations for the evolution of parton densities. As detailed numerical studies, which were 
performed in Refs. [24,26], showQ one has to take into account not only the leading pole 
terms, but also the resummed sub-leading terms, see also [28], which are not yet calculated 
completely. They are known, however, for the first two orders in a s (and partly to 0(a^)) 
for all sub- leading terms, which suggest the general form. The comparison of the leading 
and sub-leading terms shows, that at least three sub-leading terms are required to match 
the exact results. Ignorance of these contributions, as unfortunately still partly present 
in the contemporary literature, results into misleading quantitative analyzes. Sometimes 
also "DGLAP" evolution is opposed to "infrared evolution equations", etc. Here again a 
clarifying word is in order. In practice both concepts address twist-2 parton distributions. 
Their scale evolution results from the factorization of the collinear singularities and is ruled 
by the anomalous dimensions Pij(N, a s ) = YX=i a s p i 3 k 1} ( N ) ■ The corresponding Callan- 
Symanzik equations have to be solved for high enough orders in the coupling constant in the 
range of Bjorken- x, demanded by the experimental data. These equations cover the small- a; 
and the less singular terms which are equally important in quantitative analyzes. 



c For the unpolarized case see [27]. 
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3 Parton Distributions 



The polarized parton distribution functions may be determined by a QCD-analysis of the 
structure function gi(x, Q 2 ). The data analysis requires a detailed description of the denom- 
inator of the polarization asymmetry, which has to include empiric parametcrizations both 
for F 2 (x, Q 2 ) and Fl(x, Q 2 ) including potential higher twist contributions, since the region 
of Q 2 and W 2 , which is analyzed, covers rather low values, unlike the case in unpolarized 
analyzes [29]. Usually one would like to limit the data analysis to the region Q 2 <; 4GeV 2 , 
which will be possible in future measurements at a facility like EIC [30]. The present data 
sets only allow a cut Q 2 ^ lGeV 2 . In the analysis the correlation of the different pa- 
rameters of the parton distribution functions at Qq and the QCD-scale Aqcd are rather 
essential. Measuring the gluon distribution function AG(x, Q 2 ), and to some extent also the 
sea-quark distributions, the slope effects of dgi(x,Q 2 )/d\n(Q 2 ) are important. In case of 
AG there one observes a very strong correlation with a s (Q 2 ) due to the evolution equations. 
Special assumptions on Aqcd, as fixing this value to other measurements, may introduce 
severe biases. In the inclusive analysis not all the parameters chosen to model the parton 
distributions can be measured. For the data sets currently available this applies in particular 
to those parameters which describe the range of medium values of x. Their error may be- 




Figurel: The polarized parton distributions from different analyzes at the scale Q 2 = lGeV 2 , [32]. 
AAC06 : [32]; GRSV : [33]; BB : [31]; LSS : [34]. 

come rather large compared to their value. These parameters can be fixed after a first mini- 
mization and form a certain model. Their value has to be re-fitted after the global minimum 
was found, but will usually not change significantly, cf. [31] . The relative normalization of the 
different data sets is fitted within the allowed margins quoted by the experiments, to account 
for global systematic effects. In Figure 1 recent parton distribution functions [31-34] are 
compared. Further parameterizations were given in [35-38]. Within the present errors there 
is good agreement between all analyzes. The valence quark distribution functions are deter- 
mined best, with a positive polarized up-quark and and a negative down-quark distribution. 
The sea-quark distribution is found to be mainly negative, but with a larger error. To resolve 
the different flavors of the sea-quark contributions semi- inclusive data were analyzed [2] , yet 
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with rather large errors. Under certain assumptions the error on the strange-quark density 
becomes rather low [35] . Alternatively to the standard QCD fits neural- network techniques 
were used to determine the polarized parton densities in [37] . In this analysis a larger error 
than found using the conventional methods is obtained in the small x region, were data 
are sparse. If compared to earlier analyzes [31] the polarized gluon distribution comes out 
at lower values including more recent data [32,34,39]. The polarized gluon distribution 
function, although being obtained with positive central values in unconstrained fits, is still 
compatible with zero within the errors. In some analyzes [34,38,40] one demands, as second 
option, also a negative normalization of the gluon distribution in a constrained fit. The 
corresponding distribution is slightly negative and allowed by the data under the constraint 
used. The ratio AG/G was also measured in open charm photo-production [41]. Here the 
experimental errors are still large and the result is compatible with zero. Using the fit results 
of the polarized parton distributions one may form moments, cf. [31], to be compared with 
lattice simulations, in particular for the valence-quark distributions. Here the crucial point 
is to control the systematic effects and to approach realistic values of the pion mass. Cur- 
rently values in the range ~ 270 MeV can be reached in dynamical quark simulations. 
In this way new non-perturbative quantitative test of QCD will be possible soon [21,22,42]. 

4 Future Avenues 

The current picture of the polarized nucleon is still in a move and more efforts in theory 
and experiment are needed to complete it. In the forthcoming years the data-analysis 
from HERMES and COMPASS will lead to still better parton distribution functions. It 
would be important to measure the structure function F2(x, Q 2 ) at HERMES, which would 
yield an improved systematic understanding of the data. Yet the experimental precision of 
the structure function gi(x, Q 2 ) is limited and high-luminosity measurements at a future 
facility as EIC [30] is highly desirable to provide detailed QCD-tests for Aqcd, the parton 
distributions and their moments to be compared with lattice simulations. The experiments 
at RHIC will improve our knowledge on the polarized sea-quark and gluon distributions. 
Important information on the large- a; behaviour of all distribution functions, can be gained 
in the experiments at JLAB running at an increased beam energy. The HERA experiments, 
COMPASS and the JLAB experiments will finalize their measurements on deeply-virtual 
Compton scattering (DVCS) [43] during the forthcoming years and we may hope to get 
constraints on the quark angular momentum [44] using Ji's sum-rule [45]. 

As shown in [46] the theory error of the polarized gluon distribution at NLO is still large. 
The calculation of the 3-loop anomalous dimensions arc therefore required to diminish this 
uncertainty. Very high statistics measurements have to be performed in the long-term 
future for DVSC to extract constraints on the angular momentum of the gluon from the 
scaling violations of the non-forward scattering cross sections. As is well-known, the scaling 
violations of the transversity distribution function h\(x, Q 2 ) are larger than those of the non- 
singlet part of <7i (x, Q 2 ) . Detailed high statistics measurements of this quantity are desirable 
to establish this prediction of QCD experimentally Very little is known about the highcr 
twist contributions to polarized deeply-inelastic scattering. Here we may hope for results 
from the experiments at JLAB. For the general kinematic region again high-luminosity 
experiments as planned for EIC would provide an excellent opportunity. Dedicated studies 
of the twist-3 contributions to several processes should be carried out and measurements 
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shall be performed to isolate the twist-4 contributions. An interesting open issue is formed 
by twist-3 effects [4, 15] in deep-inelastic scattering in electro-weak interactions, which can 
be studied at future neutrino factories [47] . The present status of our knowledge on polarized 
parton densities is not yet sufficient and calling for refined measurements in various places 
which are crucial for the final understanding of the spin-structure of the nucleons. This will 
require extensive experimentation at a high-luminosity facility such as the future Electron- 
Ion-Collider. 
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